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In a sub-excitable light-sensitive Belousov-Zhabotinsky (BZ) chemical medium an asymmetric 
disturbance causes the formation of localized traveling wave-fragments. Under the right con- 
ditions these wave-fragment can conserve their shape and velocity vectors for extended time 
periods. The size and life span of a fragment depend on the illumination level of the medium. 
When two or more wave- fragments collide they annihilate or merge into a new wave-fragment. 
In computer simulations based on the Oregonator model we demonstrate that the outcomes 
of inter-fragment collisions can be controlled by varying the illumination level applied to the 
medium. We interpret these wave-fragments as values of Boolean variables and design collision- 
based polymorphic logical gates. The gate implements operation XNOR for low illumination, and 
it acts as nor gate for high illumination. As a NOR gate is a universal gate then we are able to 
demonstrate that a simulated light sensitive BZ medium exhibits computational universality. 
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1. Introduction 

Designing of logical gates in chemical systems can be traced back to early 1990s. Hjemfelft and col- 
leagues suggested a theoretical coupled mass flow system for implementing logic gates and finite-state 
machines|Hjelmfelt et alj [l99l] , [Hjelmfelt eralj p92] , |Hjelmfelt fc Ross[p93l , [Hjelmfelt eriill p93 



Hjelmfelt Sz Ross, 1995 and Lebender and Schneider described approaches towards building logical gates 
using a series of flow rate coupled continuous stirred tank reactors and a bistable chemical reaction [Lebender 



I&: Schneider 1994 . No experimental prototypes were implemented at that time. 

In 1994 Toth, Showalter and Steinbock presented the first ever experimental implementation of logical 
gates in the Belousov-Zhabotinsky system Toth et al, 1994 , [Toth &: Showalter 1995 . Their constructs 
of logical gates were based on the configuration of excitation wave propagation channels and the ratio 
between channel diameter and the critical nucleation radii of the excitable media. Their findings aroused 
great interest and resulted in several innovative designs of computational devices, including logical gates for 
Boolean and multiple- valued logic [Sielewiesiuk fc Gorecki 2001 , [Motoike &: Yoshikawa 2003|, Gorecki et 



2003 



al , 20091 , lYamaguchi et al , 2009 , many- input logical gates [Gorecki &; Gorecka 2006 , Gorecki &; Gorecki 



2006 , counters [Gorecki et al 2003 , coincidence detector Gorecka & Gorecki, 2003 



and distance [Gorecki et al, 2005 , Yoshikawa et al, 2009 and inductive memory 



detectors of direction 



Motoike & Yoshikawa 



All these chemical computing devices were realised in geometrically-constrained media, where exci- 
tation waves propagated along defined channels loaded with catalyst or tubes filled with BZ reagents. The 
waves perform computation by interacting at the junctions between the channels/tubes. Such an approach 
is noble, however, it essentially just imitates conventional computing architectures (wires and valves) but 
using novel materials (excitable chemical systems). Computing in unconstrained media would be a step 
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forward towards the implementation of massively-parallel chemical processors. 

To appreciate the massive-parallelism of a thin-layer chemical media we can adopt the paradigm of 
collision-based computing [Adamatzky 2003 . This paradigm originates from the computational univer- 
sality of the Game of Life [Berlekam et al , 1992] , conservative logic and the billiard-ball model [Fredkin fc 
Toffoli, 1982 and their cellular-automaton implementations [Margolus, 1984 . A collision-based computer 



employs mobile self-localized excitation to represent quanta of information in active non-linear media. In- 
formation values, e.g. truth values of logical variables, are given by either the absence or presence of the 
localizations or by other parameters such as direction or velocity. The localizations travel in space and 
collide with each other. The results of the collisions are interpreted as computation. There are no prede- 
termined stationary wires, a trajectory of the travelling localization is a momentary wire. Almost any part 
of the reactor space can be used as a wire. Localizations can collide anywhere within this space. The local- 
izations undergo transformations, form bound states, annihilate or fuse when they interact. Information 



values of localizations are transformed as a result of these collisions Adamatzky , 2003 



To implement a collision-based scheme in a spatially-extended chemical medium we must employ trav- 
elling localisations. The self-localized excitation wave-fragments, traveling in a light-sensitive BZ medium 
when it is in a sub-excitable state Sendiiia-Nadal et al 2001 are ideal candidates. These excitation wave- 



fragments behave like quasi-particles. They exhibit rich dynamics of collisions, including quasi-reflection, 
fission, fusion, and annihilation Adamatzky &; De Lacy Costello, 2007 , Toth et al', '2009] . Using the wave- 



fragments we have implemented collision-based computing schemes Adamatzky , 2004] . We have produced 



a range of basic collision-based computing schemes in computer simulations and chemical laboratory ex- 
periments including logical ga tes [Adamatzky , 2004 , Adamatzky k, De Lacy Costello , 2007 , Toth et al 



2009 , [De Lacy Costello et al, 2009 



a one-bit adder Adamatzky, 2010 



, evolvable chemical logical circuits Toth et al 2009[ , and elements of 



Adamatzky et al, 2010 



All these excitable chemical computing devices are light-sensitive, the wave-fragments grow in size 
with a decrease in illumination, and the wave-fragments collapse with an increase in illumination. In a very 
narrow illumination the wave-fragments remain localized and conserve their shape and velocity vectors 
for a (relatively) significant amount of time. Thus the wave-fragments can be used to represent quanta of 
information. What happens if the level of illumination is altered between the lower and upper limits of the 
wave-fragments-stability range? In the present paper we demonstrate that a design can be implemented 
where the outcomes of collisions between wave-fragments sensitively depend on the level of illumination. 
A logical function is realised by the collision of wave fragments. By changing the illumination level we are 



able to change the logical function. Thus we implement BZ collision-based polymorphic logical gates Stoika 
et al, 2002 , i.e. gates which change their function depending on control signals. 



The paper is structured as follows. In Sect. [2] we show how we simulate the light-sensitive Belousov- 
Zhabotinsky system. Collisions between wave- fragments are studied in Sect.|3J We describe our implemen- 
taton of collision-based polymorphic logical gates in Sect. [4} 



2. Methods 



We use the two- variable Oregonator equation Field & Noyes, 1974 adapted to a light-sensitive Belousov 



Zhabotinsky (BZ) reaction with applied illumination Beato &: Engel, 2003 



du 
dv 
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u 



{fv + 



u + q 



+ DuV^u 



u 



(1) 



The variables u and v represent the local concentrations of activator, or excitatory component, and 
inhibitor, or refractory component. Parameter e sets up a ratio of time scale for the variables u and u, q is 
a scaling parameter dependent on the rates of activation/propagation and inhibition, / is a stoichiometric 
factor. Constant (j) is the rate of inhibitor production. In the light-sensitive BZ (j) represents the rate of 
inhibitor production which is proportional to the intensity of illumination. We integrate the system using 
the Euler method with five-node Laplace operator, time step At = 0.005 and grid point spacing Ax = 0.25, 
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{a.)(t> = 0.07 



(b)(^ = 0.077 



{c)(t) = 0.07873 






{d)(t> = 0.07877 



(e)<j!> = 0.07878 



(£)<;/> = 0.079 



Fig. 1. Showing that the development of the initial excitation is sensitively dependent on the level of illumination (p. Excitation 
is initiated at the south-west edge of the vesicle. 



e = 0.022, / = 1.4, q = 0.002. The equations effectively map the space-time dynamics of excitation in the 
BZ medium and have proved to be an invaluable tool for studying the dynamics of collisions between 
travelling localized excitations in our previous work Adamatzky, 2004 , Adamatzky & De Lacy Costello 



2007 , Toth et al, 2009 De Lacy Costello et al, 2009 



The parameter (p characterizes the excitability of the simulated medium. The medium is excitable 
and exhibits 'classical' target waves, e.g. when (p = 0.07 (Fig. flk) and the medium is sub-excitable with 
propagating localizations, or wave-fragments, when is between 0.07873 and 0.07878 (Fig. flp-e). The 
medium becomes non-excitable for (/) > 0.79, and after this point wave-fragments collapse after relatively 
short time scales (Fig. [iJF) . 

When the BZ reaction is in a sub-excitable mode asymmetric perturbations lead to the formation of 
propagating localized excitation, or excitation wave-fragments. Wave-fragments of this type may travel 
in a predetermined direction for a finite period of time. If wave-fragments kept their shape indefinitely, 
we would be able to build a collision-based computing circuit of any size. In reality, the wave-fragments 
are inherently unstable: after some period of conserved-shape/ distance travelled a wave-fragment either 
collapses or expands. 

Recently [NeuNeu, 2010 , [Adamatzky et al 2010 we found a way to overcome the problem of wave- 
fragment instability via the subdivision of the computing substrate into interconnected compartments, 
so called BZ-vesicles, and allowing waves to collide only inside the compartments. Each BZ-vesicle has a 
membrane that is impassable for excitation Gorecki, 2010 , NeuNeu, 2010 . A pore, or a channel, between 



two vesicles is formed when two vesicles come into direct contact. The pore is small such that when a wave 
passes through the pore there is insufficient time for the wave to expand or collapse before interacting with 
other waves entering through adjacent pores, or sites of contact. 

A spherical compartment — BZ-vesicle — is the best natural choice as it allows for effortless arrange- 
ment of the vesicles into a regular lattice, has an almost unlimited number of input /output states and also 
loosely conforms to a structure likely to be achieved in experiments involving the encapsulation of excitable 
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chemical media in a lipid membrane |G6recki,|2010], [NeuNeu 2010|. We simulate a vesicle filled with BZ 



solution as a disc with radius R centered in (xo,yo)- Sites inside the disc are excitable, sites outside the 
disc are not excitable. We imitate wave-fragment entering the vesicle by exciting (assigning values u = 1.) 
grid nodes inside the small disc with radius r, centered in (xq + {R — s)cos{9),yo + {R — s)sin{6)). The 
following parameters are used in the illustrations: R = 100, r = 5, s = 5, 9 £ [0, 27r]. Time lapse snapshots 
provided in the paper were recorded at every 150 time steps, and grid sites with excitation level u > 0.04 
were displayed. 

3. Binary collisions 

If two wave-fragments x and y are initiated at the disc's edge at the same time they collide with each 
other, while approaching the centre of the disc. The outcome of the collision depends on the angle a 
between the velocity vectors x and y (Fig. pi). When the angle is less than some critical value /3 the 
colliding wave-fragments merge into a new wave- fragment z (Fig. ^-g) whose velocity vector is positioned 
exactly between the velocity vectors of wave- fragments x and y: z = {x + y)/2 (Fig. pk). When the angle 
between the vectors of colliding wave-fragments exceeds some critical value /3, the colliding wave fragments 
annihilate (Fig. ^-s). 

Let /3 be a critical value such that wave- fragments colliding at angle a < (3 merge into a wave- fragment, 
which propagates over an indefinitely long distance, and wave-fragments colliding at angle a > /3 annihilate. 

Proposition 1. Critical value /3 is inversely proportional to illumination level (j). 

Dependence of /3 on </> calculated in computational experiments is shown in Fig. [3^. The dependence 
is essentially linear, the deviations shown are due to digitization of the space in numerical experiments. 
Each critical value /3 has its own illumination level (j), see examples in Fig. [Sj^-f. Therefore, by varying the 
illumination of the disc we can alter the outcomes of collisions between wave fragments. 

4. Polymorphic gates 

We consider two types of gates implemented in BZ-vesicles via the collision of wave-fragments (Fig. El). 

Proposition 2. Let Boolean values of x and y be represented by wave-fragments then a BZ-vesicle imple- 
ments a two-input three-output switchable logical gate {x,y,(p) — )• {xy,x{4')^y^^y) where x(</>) = 1 (TrueJ 
if <J3 = 4>iow, O'^d (TalseJ otherwise. 

Let there be a maximum of two wave-fragments entering a BZ-vesicle. The wave-fragments enter the 
vesicle along trajectories x and y (Fig. Ek). We assume that presence of a wave-fragment at entry point 
x represents a logical value True, absence — logical value False. Similarly, if there is a wave-fragment 
entering BZ-vesicle along trajectory y we assume y=TRUE, otherwise y=FALSE. When just one of the input 
values is True then the solitary wave-fragment passes through the vesicle without significant modification 
and exits the vesicle at the site opposite its entry point (Fig. [5^-d, x = l,y = and x = 0,y = 1). If 
two wave-fragments enter the vesicle they interact and do not follow their original trajectories. Thus the 
output trajectories along which the undisturbed wave-fragments x and y move represent functions xy and 
xy, respectively (Fig. |4^). 

Interaction of wave-fragments is determined by level of illumination. When illumination is low enough, 
say 4)\ow, the colliding wave- fragments merge in a new, i.e. travelling along new trajectory, wave- fragment 
(Fig. 5i and c, x = l,y = 1). The new wave-fragment exiting the BZ-vesicle represents operation xy 
(Fig. 4j,, left). For higher level of illumination, say (phigh, the colliding wave-fragments annihilate each 
other (Fig. [5J3 and d, x = l,y = 1), no additional operation is realised. 

This is true for wave- fragments colliding at almost any angle over 7r/6. However for any particular 
angle we must select unique values of cpiow and ^high- 
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Fig. 2. Outcomes of collisions between two wave-fragments x and y approaching each other at angle a, illumination level is 
(j) = 0.07873; z shows orientation of velocity vector of the resultant wave- fragment, produced in the collision. 

Proposition 3. Let Boolean values of x, y and z be represented by wave-fragments then a BZ-vesicle imple- 
ments a three-input three-output switchable logical gate {x,z,y,(j)) — )• {xy^,xi4')z{x (B y) -\- x{^)^yz,xyz) 
where x(0) = 1 (True^ if (p = <i>iow, o^f^d (T'alsej otherwise. 



Outputs presented by trajectories of undisturbed signals x (y) — xyz [xyz) — are determined as 
follows. Wave-fragment x (y) continues traveling along its original trajectory only if neither wave-fragment 
y (x) nor wave-fragment z enter the vesicle (Fig. [4] and Fig. l6k and h, x = 1, y = 0, z = and x = 0, 
y = l,z = 0). 



6 Adamatzky, De Lacy Costello, Bull 



P 



137r/18 ^ 
IItt/IS 

Till 

7jr/18 
Su/lS 



nl6 
0.07838 0.07848 0.07858 0.07868 0.07878 0.07888 0.07898 0.07908 

(a) 




(b)« = I 



W«-S 



Fig. 3. Level of illumination determines critical angle of wave- fragments collision: (a) dependence of critical collision angle /? on 
illumination parameters 0, time lapse snapshots of wave-fragments colliding at angles /3 for their critical values of illumination 
<j>:{h)a=l,(l) = 0.0791, (c) a = ff , = 0.07883, (d) a = ^, (/> = 0.07867, (e) a = ^, = 0.07850, (f) a = ^, = 0.078. 

The following scenarios take place for both low (j)\oyj and high (j)\iigh levels of illumination. If only 
wave-fragment z is present, it travels through the vesicle undisturbed (Fig. pi x = 0, y = 0, z = 1). When 
wave-fragment z is present and also either wave-fragment x or y the wave-fragments collide and form a 
wave-fragment whose velocity vector is an average of the velocity vectors of the colliding wave- fragments. 
The newly formed wave-fragment collides with the vesicle's wall just between the output channels and 
misses both of the potential exit points. Thus no output is generated (Fig. [6],x = l,y = 0, z = l and 
X = 0,y = l,z = 1). 

For two combinations of inputs — x = l,y = l,z = and x = l,y = 1, z = 1 — the outcomes depend 
on the level of illumination. If only wave-fragments x and y or all three wave-fragments enter the vesicle 
they collide and annihilate when the level of illumination is high 0hig/i (Fig. pk). The fragments merge 
and form a new wave- fragment, which hits the output channel thus generating output value True, when 
the level of illumination is low 4>high (Fig. [6b). Thus the output channel opposite to the input channel z 
generates z{x (B y) when the level of illumination is low, and it generates xyz when the level of illumination 
is high. 

By assigning constant True to input z, we realize a two-input one-output gate {x, y, (p) — >• 
{x{<p)z{x ® y) + x{4>)^yz). Thus we arrive at the main finding of the current paper: 



Proposition 4. Fragments travelling and colliding within a BZ-vesicle implement a polymorphic logical 
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Fig. 4. Two types of logical gates controllable by illumination level 4>- 
gate switchable between functional states XNOR and NOR by changing the degree of illumination. 

5. Conclusion 

In a numerical model of the light-sensitive Belousov-Zhabotinsky (BZ) medium in a sub-excitable mode 
localized traveling excitation waves are formed. We interpreted these localizations as quanta of informa- 
tion, values of logical variables. When two or more localizations collide they annihilate or form a new 
localization. We interpreted post-collision trajectories of the localizations as the results of a computation. 
We demonstrated that by colliding wave-fragments in an encapsulated excitable chemical medium we can 
realise a number of logical gates. We showed that by changing the illumination of the chemical medium we 
could switch between different outcomes of the computation. Thus we were able to realise a polymophic 
logical gate which could execute either function XNOR or NOR depending on the level of illumination. Gate 
NOR is a universal gate, thus, as a byproduct, we demonstrated the computational universality of the BZ 
medium when in a subexcitable state. 

We would like to outline two main directions of further studies. In the theoretical part, we aim to 
focus on cascading BZ-based polymorphic gates into larger logical circuits and arithmetic schemes. In 
experimental part, we aspire to implement theoretical constructs in chemical laboratory experiments. 
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Fig. 5. Implementation of polymorphic gate {x,y,(j>) — )■ {xy , x(4>)xy ,xy) . Scheme of the gate is shown in Fig. Bk. Time 
lapse snapshots of the sub-excitable media are shown for various illumination levels <j) and collision angles ct. Wave-fragments 
represent logical values of inputs x and y. Inputs (entry points, pores) are marked by thin lines, outputs (exit points, pores) 
are marked by thick lines. 
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Fig. 6. Implementation of polymorphic gate {x, z, y, 0) — >■ {xyz, xi<P)z{x (By) + xW^V^t xyz). Scheme of the gate is shown 
in Fig. Wn. Time lapse snapshots of sub-excitable medium are shown for various illumination levels (f> and collision angles a. 
Waves represent logical values of inputs x and y. Wave-fragments x and z, and z and y collide at angle ^; wave-fragments x 
and y collide at angle 5. Inputs (entry points, pores) are marked by thin lines, outputs (exit points, pores) are marked by 
thick lines. 
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